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NOL  HYPERSONIC  TUNNEL  NO.  4 RESULTS  V. 
EXPERIMENTAL  AND  THEORETICAL  INVESTIGATION  OF 
A COOLED  HYPERSONIC  WEDGE  NOZZLE 

By 

P.  Wegener,  R.  K.  Lobb,  E.  M.  Winkler, 

M.  Sibulkin  and  H.  Staab 


ABSTRACT:  A new  w?ter  cooled,  wedge-type  nozzle  has 

been  in  operation  in  the  NOL  12  x 12  ca  Hypersonic 
Tunnel  No.  4 since  January  1952.  This  report  presents 
results  of  theoretical  and  experimental  investigations 
covering  the  performance  of  this  nozzle.  It  was  found 
from  measurements  of  pitot  and  static  pressure  that  the 
nozzle  produced  a shock-free,  almost  isentropic  expansion 
having  a test-section  Mach  number  distribution  comparable 
in  quality  to  that  obtained  in  existing  conventional 
nozzles  at  lower  liach  numbers.  The  cooling  system  used 
was  adequate  to  maintain  the  nozzle  block  temperatures 
at  conotant  low  values.  Tha  temperature  distribution 
in  a nozaple  block  was  also  measured  and  it  showed  a 
rise  near  the  throat.  This  rise  agreed  qualitatively 
with  theoretical  calculations.  A comparison  of  the 
calculated  and  measured  rate  of  heat  transfer  to  a 
nozzle  block  suggested  that  boundary  layer  transition 
occurs  after  the  throat.  A pressure  and  temperature 
survey  at  the  end  of  a nozzle  block  where  the  boundary 
layer  is  turbulent  showed  that  the  velocity  distribution 
is  similar  to  that  found  at  low  speed.  This  survey  was 
also  used  to  calculate  the  overall  heat  transfer  rate  to 
the  nozzle  block.  This  result  agreed  with  that  obtained 
from  the  direct  method  of  determining  the  heat  transfer 
by  measuring  the  cooling  water  rate  and  temperature 
increase . 

In  conclusion  it  is  found  that  a hypersonic  wedge  nozzle 
produces  shock  free  expansions  and  acceptable  flow  in  the 
test  section  and  serves  well  as  a research  tool. 


U.  S.  NAVAL  ORDNANCE  LABORATORY 
WHITE  OAK,  MARYLAND 
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carried  out  in  the  continuous  NOL  12  x 12  o Hypersonic 
Tunnel  No.  4.  This  facility  was  first  put  into  operation 
in  May  1950.  The  titles  of  the  previous  NAVORDS  discussing 
results  froa  the  tunnel  are: 

I Air  Liquefaction.  NAVORD  1742,  4 Jan  1951 

II  Diffuser  Investigation.  NAVORD  2376,  5 May  1952 

III  Diffuser  Investigation  with  Models  and  Supports, 
NAVORD  2435,  1 July  1952 

IV  High  Supply  Temperature  Measurement  and  Control, 
NAVORD  2574,  8 Oct  1952 

The  present  NAVORD  (Results  V)  presents  an  experimental 
and  theoretical  discussion  of  a cooled  high  Mach  number 
noxzle.  The. work  was  jointly  sponsored  by  the  U.  S.  Naval 
Bureau  of  Ordnance  and  the  U.  S.  Air  Force,  Flight  Research 
Laboratory. 

The  authors  are  indebted  to  Messrs.  Charles  B.  rhite  and 
E.  J.  Stollenwerk  for  the  mechanical  design  of  the  nozzle. 
Prof.  n.  naxiJs  advice  on  the  cooling  problem  is  gratefully 
acknowledged , Messrs.  L-.  Liccini  and  R . Garren  participated 
during  the  tests  and  Messrs.  J.  Kendall,  Jr.  and  H.  Maxwell 
assisted  in  making  the  boundary  layer  calculations. 


EDVARD  L.  WOODYARD 
Captain,  U3N 
Commander 


H.  H.  KURZVEG,  Chief 
Aeroballistic  Research  Department 
By  direction 
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I ILLUSTRATIONS 

Figure  1.  The  NOL  12  x 12  ca  Hypersonic  Tunnel  No.  4. 

' Figure  2.  Cooled  Wedge  Nozzle  Block  with  Thermocouple 

Plugs  Rescued  * 

Figure  3.  Cooled  Wedge  Nozzle  ior  Tunnel  No.  4. 

Figure  4.  Scblieren  Photograph  of  the  Flow  Near  the 
Throat  o£  a Minimum  Length,  M m 5.18, 

Nozzle,  pQ  s 4.5  atm,  TQ  = 3330Ko 

Figure  5.  Shadowgraph  and  Schlieren  Photographs  of 
Flow  lu  the  Wedgo  Nozzle . 

Figure  6.  Mach  Number  Variation  in  Nozzle  Exit  vs 

Time  for  Changing  Coolant  Hate  pQ  x 10  atm, 

T - 593°K. 

Figure  7.  Static  Pressure  Distribution  along  Centerline 
of  Nozzle  Sidewall,  M - 7.0,  pQ  ■ 10  atm, 

T0  s 593°K. 

Figure  8.  Mach  Number  Distributions  along  Surveys  in  a 
Plane  One  Inch  Upstream  from  the  Nozzle  Exit, 
pQ  5 10  atm,  Tfl  s 593°K. 

Figure  9.  Comparison  of  measured  and  Calculated  Boundary- 
Layer  Growth  along  Nozzle  Wall. 

! Figure  10.  Measured  Turbulent  Boundary-Layer  Profiles  with 

Heat  Transfer  at  a Free  Stream  Mach  Number  of  7. 

Figure  11 . The  u + , y * Representation  of  the  Measured 
Velocity  Profiles  at  S s 7 and  Comparison 
with  Theory. 

Figure  12.  Temperature  Distribution  in  Nozzle  InlGt  at 
Three  Different  Levels  of  Supply  Temperature 
(Temperature  in  °K). 

i Figure  13.  Nozzle  Wall-Temperatures  Near  Throat  (Station 

No.  13)  at  Various  Supply  Temperatures , 

‘ i pa  a 10  atm,  M = 7.6.  Cooling  Rate:  6.5 

» Gallons /Min.  * 

! Figure  14.  Nozzle  Wall-Temperatures  Near  Throat  (Station 

I No.  13)  at  Various  Supply  Pressures,  T^  z 610°K, 

Mg  s 7.6.  Cooling  Bate:  6.5  Gallons/Kin. 

i Figure  15.  Wall  Temperatures  at  9 Stations,  Temperature 

| Gradient  in  Wall,  and  Extrapolated  Surface 

! Temperatures  of  Nozzle  ML  ■ 7.6,  p0  s 10  atm, 

1 Tc  ■ 593°K.  Cooling  Rate:  6.5  Gallons/Mic. 

i Figure  16.  Wall  Temperatures  at  S Stations,  Temperature 

\ Gradient  in  Wall,  and  Extrapolated  Surface 

i f Temperatures  of  Nozzle,  M_  z 7.6,  pp  = 20  ata, 

| \ T0  s 593°K.  Cooling  Rate:  6.5  Gallon* /kin. 

it  Figure  17.  Isotherms  in  Nozzle  for  Hj  s 7.6,  Tn  a 593°K, 

if  p0  i 1C  and  20  atm,  Cooling  Rate:  6.5  Gallons /Min . 

i V 


iv 

UNCLASSIFIED 


1 


UNCLASSIFIED 
NAVORD  Report  2701 


Figure  IS . 


Figure  19. 


Figure  20. 


Figure  21. 


Figure  22. 
Figure  23. 
Figure  24. 
Figure  25. 
Figure  26. 
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Figure  28. 


Figure  29. 


Figure  30. 


Figure  31. 


Figure  32. 


Figure  33. 


Figure  34. 


Nozzle  Wall-Temperatures  as  Function  of 
Running  Tine  at  Boundary-Layer  Survey  Station, 

Kg  * 7.6,  p-  .=  21.4  atm,  T = 593°R.  Cooling 
Rate:  6.5  Gallons/Min. 

Temperature  Drop  in  n n**l#  Wall,  = 7.6, 
pQ  = 21.4  fata,  TQ  = 593°K.  Cooling  Rate: 

6.5  Gallons/Min, 

Temperature  Increase  of  Cooling  Wafer  as 
Function  of  the  Cooling  Rate  for  Tunnel 
Operation,  = 8.0,  p0  s 20  atm,  TQ  r 600°K. 

Heat  Flo*  Rate  to  Cooling  Water  as  Function 
of  Cooling  Rate  for  Tunnel  Operation, 

M = 8.0,  ?0  r 20  ata,  T = 600°K. 

Coolant  Tube  Position?  and  Distensions  Used  in 
Heat-Transfer  Analysis. 

Variation  of  k/ke  with  Coolant  Tube  Size  and 
Position. 

Deteraining  Nozzle  Surface  and  Coolant  Tube 
Surface  Temperatures . 

Velocity  Distribution  along  Nozzle,  K = 8.0, 

T0  = 600°K.  8 

Acceleration  Distribution  along  Nozzle, 

Mg  a 8.0,  T = 600°S. 

Heat-Transfer  Coefficient  Distribution 

Assuming  Laminar  Flo*  and  Neglecting 

the  Effects  of  Pressure  Gradients,  Mg  • 5.0, 

T„  = 600°K,  p„  « 20  atm. 

ESat-Transfer''Coefficient  Distribution 
Assuming  Turbulent  Flo*  and  Neglecting  the 
Effects  of  Pressure  Gradients,  = 8.0, 

Tq  - 600°K,  p0  - 20  at*. 

Growth  of  Momentum  Thickness  Assuming  Laminar 
Flo*  and  Including  the  Effects  of  Pressure 
Gradients,  Hg  ■ 8.0,  T s 600°S,  p,,  = 20  ata. 
Heat-Transfer  Coefficient  Distribution  Assuming 
Laminar  Flo*  and  Including  the  Effects  of 
Pressure  Gradients,  Mg  = 8.0,  T0  - 600°K, 
p s 20  at*. 

Growth  of  Momentum  Thickness  Assuming  Turbulent 
Flo*  and  Including  the  Effects  of  Pressure 
Gradients,  Mg  ; 8.0,  T : 600°K,  p m 20  atm. 
Heat-Transfer  Coefficient  Distribution  Assuming 
Turbulent  Flo*  sad  Including  the  Effects  of 
Pressure  Gradients,  Mg  - 8.0,  T - 600°K,  pG  a 20  ata. 
Comparison  of  Heat-Transfer  Coefficients  Showing 
Effect  of  Type  of  Boanda.ry-Lp.yer  and  Pressure 
Gradients.  Mg  x 8.0,  T = 600^,  p„  s 20  stss. 
Installation  of  Thermocouple  Plug  in  Nozzle  Block. 
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NOL  HYPEB30NIC  TUNNEL  NO.  4 RESULTS  V: 
EXPERIMENTAL  AND  THEORETICAL  INVESTIGATION  OF 
A COOLED  HYPERSONIC  WEDGE  NOZZLE 


INTRODUCTION 


1.  During  the  initial  planning  phase  (reference  1} 

of  the  NOL  12  x 12  ca  Hypersonic  Wind  Tunnel  (Figure  1), 
it  was  decided  to  use  two  plane  blocks  as  a nozzle 
(Figures  2 and  3)  producing  a continuous  expansion  of  the 
air.  The  Mach  number  in  any  section  of  such  a "vedge 
nozzle",  and  in  particular,  the  maxiaun  Each  nusber  in 
the  exit  section  can  then  be  easily  varied  by  varying 
the  nozzle  area  ratio  only.  This  feature  appeared  to  be 
advantageous  for  the  study  of  such  basic  phenomena  as 
aii  condensation  (reference  -i) , needed  pressure  ratios 
to  operate  at  high  Mach  numbers  (reference  39),  etc. 

It  was  planned  to  later  attempt  the  design  of  a nozzle 
producing  uniform  flow  in  a greater  region  fur  testing 
larger  models.  Another  reason  for  the  initial  selection 
of  such  a simple  nozzle  was  the  fact  that  boundary  layer 
corrections  of  theoretical  potential  flew  nozzle  walls 
designed  to  produce  such  uniform  floors  are  insufficiently 
known  at  high  2ach  numbers. 

2.  Experimentation  with  t**o  types  of  wedge  nozzles 
during  the  last  2-1/2  years  has  shown  that  these  nozzles 
produce  satisfactory  flow,  and  a theoretical  and  experi- 
mental discussion  of  the  nozzle  performance  with  some 
comparison  with  uniform  flow  nozzles  will  be  given  in 
section  XI  of  this  report. 

3.  To  prevent  air  condensation  in  the  wind  tunnel  at 
high  Mach  numbers,  the  supply  air  must  be  preheated.  Even 
highly  purified  air  does  not  supercool  appreciably  (ref- 
erence 2)  in  contrast  to  the  well  known  case  of  rapid 
expansions  of  water  vapor  in  nozzles  (reference  3) . It 

is  therefore  advisable  that  the  expansion  of  the  air  in 
the  tunnel  does  not  enter  the  condensation  region.  Wind 
tunnel  supply  temperatures  and  pressures  for  such  operation 
are  given  on  Figure  8 of  reference  4.  Although  these 
temperatures  are  much  lower  than  those  which  would  simulate 
actual  free-flight  stagnation  temperatures  (compare  Figure  1, 
reference  5) , they  are  high  enough  to  pose  a severe 
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technical  problem  at  Each  nuabers  above  say  E - 8. 

At  the  peak  Each  nuabers  for  practicable  tunnel  oper- 
ation, one  aay  therefore  well  resort  to  extreme 
purification  to  obtain  scss  super-cooling  and  resultant 
saving  in  heating  as  pointed  out  by  Nagaaatsu  (referonce  2) . 
In  this  report  ve  Till  deal,  however,  aainly  with  Mach 
nuabers  of  the  order  7 to  8 where  a moderate  supply 
pressure  range  froa  5 to  30  atmospheres  at  temperatures 
of  the  order  300  to  400°C  sakes  continuous  tunnel 
operation  quite  practical. 

4.  During  operation  at  elevated  supply  temperatures, 
all  components  of  the  tunnel  are  being  heated  by 
conduction,  forced  convection  and  radiation.  The  nozzle 
and  test  section  walls,  for  example,  will  tend  to  reach 
a temperature  close  to  the  recovery  temperature  after 
some  time  of  operation.  This  recovery  temperature  is 
in  turn  of  the  order  of  the  temperature.  Depending 

on  the  local  rate  of  heat  transfer  from  the  flow,  the 
locally  different  heat  capacity  of  the  tunnel  components, 
etc . , the  equipment  will  become  non-unif ormly  hot . An 
example  of  such  uneven  heating  of  the  nozzle  surface  and 
other  components  is  given  on  Figures  4 and  6 of  reference  5. 
It  can  be  seen  that  the  throat  of  the  nozzle  heats  up  much 
mere  rapidly  than  the  nozzle  end  (compare  section  V).  If 
supply  p>  assure  and  temperature  are  held  constant  for  a 
fixed  Each  number,  one  expects  that  finally  all  tunnel 
components  will  reach  some  constant  temperature.  Such  a 
constant  temperature  could  not  be  attained  in  our  tunnel 
after  several  hours  of  operation.  On  the  other  hand,  the 
slowly  Increasing  temperatures  of  all  tunnel  components 
cause  warpage  of  metal  parts,  breakdown  of  all  rubber 
based  seals,  glass  breakage,  etc.  Warpage  is  particularly 
effective  at  the  narrow  nozzle  throat  where  small  area 
changes  result  in  major  Each  number  changes  in  the  test 
section,  (see  section  II).  In  order  to  avoid  these 
difficulties,  retain  the  flexibility  of  Each  number  change 
inherent  in  a two-dimensional  wedge  nozzle,  and  rapidily 
attain  a constant  temperature  and  geometry,  a nozzle 
water-cooling  system  was  designed  and  built,  (Figures  1, 

2 and  3).  It  was  further  found  that  the  top  of  the  test 
section  frame  had  to  be  cooled  (Figure  1)  to  obtain  fixed 
base  points  for  the  nozzle  bed.  With  such  an  arrangement, 
the  following  is  found  (see  sections  XZ  and  IV): 
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(a)  Most  of  the  nozzle  surface  remains  near 
room  temperature  at  all  times.  (In 

our  tunnel  tho  cooling  water  discharge  is 
of  the  order  10  gallons  per  minute.) 

(b)  Constant  Mach  number  readings  in  the 
test  section  are  attained  after  a few 
minutes  of  operation. 

(c)  The  tunnel  can  be  operated  continuously 
for  hours. 

(d)  Materials,  gaskets,  etc.  do  not  undergo 
permanent  changes  and  ail  operating 
conditions  can  be  repeated  accurately. 

5.  The  following  sections  of  this  NAVORD  report  give 
an  account  of  some  of  the  theoretical  and  experimental 
work  done  to  date  on  this  nozzle. 


II.  NOZZLE  DESIGN  AND  PERFORMANCE 

6.  Design  Problems  of  Hypersonic  Nozzles.  The  aero- 
dynamic design  of  a potential  flow  nozzle  for  Mach 
numbers  above  5.0  is  the  same  as  that  for  conventional 
supersonic  nozzles.  There  are  several  methods  (all 
based  on  theory  of  characteristics)  for  determining 
the  contours  of  both  two  and  three-dimensional  nozzles. 

A two-dimensional  nozzle  can  be  designed  bv  the  graphical 
method  of  Frandtl  and  Busemam  (reference  6)  or  by  the 
semi-graphical  method  given  by  Puckett  (reference  7 and  8). 
These  methods  are  tedious  and  are  subject  to  graphical 
error  especially  in  the  hypersonic  range  for  which  the 
characteristic  net  becomes  greatly  elongated.  The 
analytical  methods  of  Foelsch  (reference  9)  and  Atkin 
(reference  10)  are  usually  preferred  for  practical  reasons 
because  they  are  more  accurate  and  can  be  used  on  high 
speed  computing  machines.  In  the  three-dimensional  case 
the  design  is  complicated  by  the  fact  that  the  characteristic 
net  in  the  hodograph  plane  has  to  be  determined  step  by 
step  aloag  with  the  net  in  the  physical  plane.  However, 
for  axially  symmetric  flows  there  are  available  several 
graphical,  numerical  and  analytical  methods  (see  references 
11,  12  and  13).  A three-dimensional  nozzle  of  arbitrary 
cross-section  can  be  obtained  by  tracing  the  stream  lines 
that  outline  the  area  throughout  a known  axially  symmetric 
flow  (see  references  14  and  15) . This  method  has  the 
disadvantage  that  the  cross-sectional  shape  so  obtained 
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varies  along  the  nozzle  thus  making  it  difficult  to 
fabricate.  It  may  be  possible  to  correct  the  surfaces 
to  give  a constant  cross-sectional  shape  by  a linearized 
method  of  characteristics  such  as  given  by  Ferri  in 
reference  16.  Another  type  of  three-dimensional  nozzle 

AMnl  AcH  ncr  o cniiorrt  f hi*AO  ♦ anH  ovnen/tlmy  -i  n r. 

r - “ o — 1.  — ——  — - - - — — — v**w  W VK.^wty  ) 

first  in  one  plane  then  in  the  plane  at  right  angles, 
is  described  in  reference  17. 

7.  one  of  the  main  practical  problems  in  hypersonic 
nozz?e  design  is  the  achievement  of  a smooth  flow  through 
the  nozzle  throat.  For  high  Mach  numbers  the  ratio  of 
the  test  section  area  to  throat  area  becomes  very  large . 

As  a result  for  medium  sized  tunnels  the  throat  must  be 
a narrow  slit  or  some  configuration  of  equivalent  area. 
This  is  illustrated  in  the  table  below.  The  throat 
must  therefore  be  machined  very  accurately.  The 
mechanical  design  is  further  complicated  by  the  possible 
warpage  of  the  throat  due  to  the  high  surface  temperatures 


1 


m 


TABLE  I I -I 

Throat  and  Exit  Areas  for  Two  and  Three-Dimensional 

(M  = 8)  Nozzles. 


Two-dimensional 


///  //  12  x 12  cm 


T'hf  AO*- 


.063  x 12  cm 


' T 

-/i| 


.Ot— . 


•„vL 


Axially  Symmetric 


Square  throat  and 
exit  plane 


12  cm  dia.  0 .87  cm  dia. 


12  x 12  cm 


.87  x .87  cm 


that  can  occur  in  this  region  during  "hot”  tunnel  operation 
(see  section  IV  and  V).  Also  the  subsonic  section  isust 
usually  withstand  high  pressures  while  in  the  supersonic 
portion  the  pressure  is  very  low.  (The  NOL  hypersonic 
tunnel  nc.  4 operates  at  supply  pressures  up  to  30  ate). 
High  supply  pressures  are  required  in  order  to  have 
continuum  flow  in  the  test  section  at  high  Mach  numbers. 
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8.  A two-dimensional  nozzle  was  selected  for  the  NOL 
hypersonic  tunnel  for  the  following  reasons.  (1)  Nozzle 
boundary-layer  corrections  are  r.ot  accurate  at  high  Mach 
numbers  (primarily  due  to  the  lack  of  skin-friction  data, 
see  reference  15).  In  th«  t.wo-dimes5  icnal  case  the 
boundary  layer  can  be  allowed  for  approximately,  by  ad- 
justing the  nozzle  blocks.  For  an  axially  symmetric  or 
square  nozzle,  on  the  other  hand,  errors  in  the  boundary- 
layer  correction  resulting  in  poor  Mach  number  distribution, 
would  require  reaachining  the  surfaces  or  discarding  the 
nozzle.  (2)  Disturbances  in  an  axially  symmetric  nozzle 
tend  to  focus  along  the  axis.  (3)  A two-step,  three- 
dimensional  nozzle  is  unsatisfactory  (see  reference  1?) 
because  of  boundary-layer  affects.  (4)  The  cost  of 
fabricating  a three-dimensional  nozzle  is  considerably 
higher  than  for  a two-dimensional  nozzle. 


9.  For  the  initial  tests  in  the  hypersonic  tunnel  a 
two-dimensional  wedge  nozzle  (straight  diverging  walls) 
was  chosen  in  preference  to  a uniform  flow  nozzle. 

With  this  wedge  nozzle  a wide  range  of  Mach  numbers  can 
be  achieved  by  simply  changing  the  throat  opening  without 
introducing  disturbances  due  to  incorrect  nozzle  contour. 
The  flow  in  a hypersonic  wedge  nozzle,  for  small  wall 
angles,  closely  approximates  two-dimensional source  flow 
(constant  properties  along  circular  arcs  with  centers 
roughly  at  the  throat).  In  a uni for™  flow  nozzle,  on  tbs 
other  hand,  the  flow  expands  slowly  along  the  nozzle 
walls  as  compared  to  that  along  the  axis.  This  produces 
relatively  strong  "cross-flow"  pressure  gradients  which 
cause  the  boundary  layer  on  the  sidewalls  to  build  up 
along  the  centerline  (references  19  and  20).  Also  a 
wedge  nozzle  of  the  same  length  and  Mach  number,  as  a 
uniform  flow  nozzle,  has  a smaller  maximum  wall  angle. 
This  reduces  the  possibility  of  flow  separation. 


10.  Previous  investigations  in  a minimum  length  uniform 
flow  nozzle  designed  for  M = 5.18  have  shown  that  a small 
radius  at  the  throat  can  lead  to  such  flow  separation  and 
shock  waves  (see  Figure  4).  Here  the  throat  radius  is 
about  1/16"  and  the  wall  angle  slightly  downstream  of  the 
throat  is  36.2°.  In  order  to  minimize  this  effect,  the 
throat  of  the  wedge  nozzle  was  formed  as  a cylinder  of  3/8” 
radius  and  the  wall  angle  was  kept  to  approximately  5-1/2°. 
Shadowgraph  and  schlieren  photographs  of  the  flow  near  the 
throat  of  the  wedge  nozzle  for  a Mach  number  setting  of  5.0 
are  given  in  Figure  5.  (This  setting  corresponds  to  the 
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largest  throat  opening  possible  with  the  present  nozzle). 
From  these  photographs  we  can  see  that  there  is  no 
separation  oe  large  disturbances  in  the  flow  field.  The 
triangular  pattern  appearing  in  the  nchlieren  photograph 
agrees  rovghly  with  that  found  graphically  by  the  method 
of  characteristics.  Tbs  distortion  of  the  throat  in  the 
shadowgram  is  optical  and  is  due  to  the  large  density 
gradients  in  this  region  of  flow. 

11.  When  the  tunnel  is  operated  at  high  supply  temperatures 
the  nozzle  blocks  become  heated  and  tend  to  expand,  par- 
ticularly in  the  throat  region  where  the  rate  of  heat 
transfer  is  large.  The  expansion  is  inward  because  the 
tunnel  casing,  to  which  the  blocks  are  attached,  remains 
relatively  cold.  As  a result  the  throat  area  decreases 

and  the  nozzle  exit  Mach  number  increases  with  time.  By 
introducing  a cooling  system  (see  Figures  1,  2 and  3)  the 
nozzle  dimensions  can  be  kept  constant  after  a short  period 
of  operation  (For  p0  = 10  atm,  T = 593QK  and  M a 7.0  the 
nozzle  blocks  reach  equilibrium  after  about  ten  minutes  of 
operation) . The  sensi tivity  of  Mach  number  to  coolant  rate 
is  illustrated  in  Figure  C.  It  can  be  seen  that  without 
any  nozzle  cooling  the  Mach  number  increases  with  time. 
Nozzle  cooling  also  -tends  to  keep  the  boundary-layer 
characteristics  constant  with  time. 

12.  Flow  Calibration.  An  example  of  the  flow  calibration 
will  Kft  Vlpnn  f<vr  a Barh  mimhor  n£  7.0  CSl" . HoWSVCr , the 
results  are  generally  true  for  other  Mach  number  settings 
This  particular  setting  corresponds  to  that  chosen  for  the 
boundary-layer  survey  and  the  heat-transfer  measurements 
given  in  sections  III  and  XV. 

13.  The  flow  in  a plane  1"  upstream  from  the  nozzle  exit 
was  surveyed  with  pitot  and  static  pressure  probes.  Static 
pressures  along  the  centerline  of  the  nozzle  siderwall  were 
also  measured.  The  pitot  probes  were  made  from  .095" 
hypodermic  tubing  with  a wall  thickness  of  .01".  The  static 
pressure  probes  were  made  from  the  same  tubing  with  an  8° 
cone  attached  to  the  tip.  Four  .025"  orifices  were  located 
20  probe  diameters  aft  of  the  shoulder.  The  static  orifices 
in  the  nozzle  sidewall  hcd  an  opening  of  025".  In  order  • 
to  measure  the  low  static  pressures  near  the  nozzle  exit 

(of  the  order  of  2ma  Hg)  it  was  necessary  to  use  a silicone 
oil  micro-manometer  equipped  with  an  optical  reading  system. 
This  instrument  has  an  accuracy  of  better  than  1/2%  in  the 
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range  under  consideration.  A precision  mercury  manometer 
or  calibrated  aneroid  type  gauges  were  used  for  measuring 
the  higher  pitot  pressures. 

14.  A plot  of  the  static  pressure  distribution  along 
the  centerline  of  the  nozzle  sidewall  in  terms  of  the 
supply  pressure,  for  an  exit  Mach  number  of  7.0,  is  given 
in  Figure  7.  There  appears  to  be  a smooth  expansion  along 
the  entire  length  of  the  nozzle.  The  pressure  gradient  at 
the  nozzle  exit  corresponds  to  a drop  in  Mach  number  of 
about  .02  M/cm.  This  compares  with  the  value  of  .03  M/cm 
derived  from  one-dimensional  theory  assuming  inviscid 
flow  throughout  the  nozzle  in  which  case  the  growing 
boundary  layer  on  the  nozzle  wall  is  neglected. 

15.  The  Mach  number  distribution  for  various  traverses 
in  a plane  1"  upstream  from  the  nozzle  exit  as  determined 
from  pitot  and  static  pressure  measurements  at  the  identical 
location  using  the  Rayleigh  formula  are  given  in  Figure  8. 
The  maximum  variation  in  Mach  number,  outside  of  the 
boundary  layer,  is  2%.  In  earlier  surveys  it  was  found  that 
there  were  disturbances  in  the  flow  origination  near  the 
junctions  of  the  nozzle  throat  with  the  sidewalls.  By 
using  metal  to  metal  g'-als  between  the  nozzle  and  sidewalls 
these  disturbances  could  be  practically  eliminated  as 
indicated  by  the  smooth  Mach  number  distributions.  From 
these  measurements  we  find  the  stagnation  pressure  loss 
throughout  the  nozzle  to  be  about  4%.  A schlieren 
photograph  of  the  flow  in  the  nozzle  exit  for  a Mach  number 
of  7.0  is  given  in  Figure  5. 

16.  Boundary-Layer  Growth . If  at  some  station  along  the 
nozzle  the  Mach  number  in  the  free  stream  is  known  then 
ve  can  estimate  the  boundary-layer  displacement  thickness 
on  the  nozzle  walls  by  comparing  the  flow  area  obtained 
from  isontropic  flow  tables  with  the  geometrical  area 
enclosed  by  the  walls.  (A  one-dimensional  analysis  is 
sufficient  here  since  it  closely  approximates  the  two- 
dimensional  analysis  for  a wall  angle  of  only  3-1/2°.) 

In  the  calculations  we  assume  that  the  boundary  layer  is 
of  equal  thickness  on  all  four  walls  and  zero  at  the 
throat  (In  section  Y the  displacement  thickness  is  sho^n 
to  be  about  . o05  mm  at  the  throat).  Also  since  the 
stagnation  pressure  loss  in  this  nozzle  is  small,  we  can 
determine  the  flow  area  with  good  accuracy  from  simply 
the  ratio  of  the  static  pressure  to  the  supply  pressure 
using  iuentropic  flow  tables.  Therefore,  for  each  static 
pressure  measurement  on  the  nozzle  sidewall  we  can  obtain 

7 

UNCLASSIFIED 


UNCLASSIFIED 
NAVORD  Report  2701 


a displacement  thickness.  The  growth  of  tbe  boundary 
layer  displacement  thickness  along  the  nozzle  obtained 
in  this  manner  is  compared  to  that  calculated  by  Tucker's 
method  (reference  18)  in  Figure  9.  In  the  theory  the 
walls  are  assumed  to  be  thermally  insulated,  whereas  the 
actual  nozzle  walls  are  cooled.  This  probably  accounts 
for  the  experimental  values  being  thinner  (i.e.  increased 
density  near  the  wall)  than  those  calculated  by  the 
theory.  The  displacement  thickness  determined  directly 
from  the  boundary- layer  survey  (see  section  III)  fits 
well  with  the  other  experimental  points.  If  is  interesting 
to  note  that  the  total  thickness  of  the  boundary  layer  at 
the  nozzle  exit  obtained  by  Tucker's  method  is  24  mm 
which  is  close  to  the  26  wja  measured  in  the  boundary  layer 
survey . 


III.  TURBULENT  BOUNDARY-LAYER  CHARACTERISTICS 
AT  END  OF  NOZZLE 

17.  In  addition  to  the  indirect  determination  of 
growth  of  the  displacement  thickness  along  the  nozzle 
wall  discussed  in  the  last  section,  a detailed  survey 
of  the  boundary-layer  profiles  at  one  point  was  undertaken. 

A survey  was  made  perpendicular  to  the  west  nozzle  wall  one 
inch  from  the  nozzle  exit.  The  data  were  unaffected  by  the 
juncture  of  nozzle  end  and  first  diffuser  plates.  Since 
the  position  of  the  transition  point  along  the  nozzle  is  still 
unknown,  it  is  interesting  to  note  that  at  this  survey  station 
a fully  developed  turbulent  boundary  layer  was  found.  The 
pertinent  results  are  summarized  in  the  following  table. 


TABLE  III:  Boundary-Layer  Survey 

Station:  Center  of  west  nozzle  wall,  1”  upstream  of 

' nozzle  exit. 

Supply  pressure:  ■ 10  atmospheres 

Supply  temperature:  T0  = 593°K  (50°  above  that  temperature 

needed  to  avoid  air  condensation 
throughout . ) 

Wall  temperatui  a at  survey  station:  Tw  = 316°k  (cooling 

system  in  operation) 

Free  stream  temperature:  T«  = 56°K  (Tw/T«©  = 5.7) 

Free  stream  Mach  number:  lLo  = 7.0  t 0.01 
Free  stream  velocity:  = 1040  m/sec 

Free  stream  Mach  number  gradient:  appr.0.02  M/cm  (compare 

section  II) 
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Reynolds  nuabei*  based  on  displacenent  thickness: 

Re  (S*  ) s 43,000 

Re  based  on  equivalent  flat  plate  length  (calculated) : 

Re  (x)  « 107 

Boundary-layer  thickness  at  point  where  XK.  -a  0.995\A<g» , 

| ; 26  n 

Displacement  thickness = 10.5  me 
Momentum  thickness:  & : 1.24  n 

Energy  thickness  (see  section  ?):  <p  = 2.78  m®. 

18.  The  above  data  were  taken  in  the  following  manner: 
Static  pressure  was  measured  by  a hole  1/2  sm  diameter 
in  the  wall  and  by  a 5°  cone-cylinder  static  probe  in 
the  free  stream  just  outside  the  edge  of  the  boundary 
layer.  (Silicone  oil  manometers  of  about  i 2 microns 
measuring  accuracy  during  tunnel  operation  were  used.) 

Since  agreement  between  these  two  static  pressure  was 
found  to  be  of  the  order  of  1 per  cent,  a constant  value 
of  static  pressure  was  assumed  to  exist  through  the 
boundary  layer. 

19.  Pitot  pressure  was  surveyed  from  wall  to  free  stream 
with  a flattened  hypodermic  tube  of  about  0.18  mm  opening. 
This  opening  is  large  enough  to  avoid  errors  due  to  slip 
flow  effects  on  the  impact  pressure  measurement  in  the 
region  of  low  Reynolds  numbers  in  the  boundary  layer  near 
the  wall  (reference  21). 

20.  From  measured  pitot  pressure  and  static  pressure, 
the  local  Mach  number  could  be  determined  from  the 
Rayleigh  formula  (reference  8 p.  95).  The  Rayleigh  Mach 
number  is  shown  on  Figure  10a.  (The  assumption  of  the 
validity  of  the  perfect  gas  law  implicit  in  this  evaluation 
and  the  use  of  K • 1.4  is  warranted  in  our  range  of 
thermodynamic  data  (reference  22,  Figure  19).  The  distance 
from  the  waliy , is  made  dimensionless  by  the  boundary- 
layer  thickness  § « 26  mm.  At  the  measured  point  closest 
to  the  wall  (Vs  mm),  the  flow  is  still  supersonic 

( M * 1.3).  * 


21.  Total  temperature  was  surveyed  with  a shielded, 
total  temperature  probe.  This  probe  was  calibrated  in 
the  free  stream  in  approximately  the  range  of  Mach  numbers, 
densities  and  temperatures  encountered  in  the  boundary- 
layer  measurements.  The  probe  recovery  factor 
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is  smaller  than  one  because  of  instrumental  errors  due 
to  radiation  and  conduction  from  the  thermocouple  junction 
to  the  surroundings.  Such  errors  are  particularly  hard 
to  overcome  at  the  high  Mach  numbers  (lorn  test  section 
density)  and  high  temperatures  (large  radiation  losses) 
in  question.  Using  the  measured  free  stream  recovery 
factor  to  correct  the  boundary-layer  total  temperature 
measurement,  ve  obtain  Figure  10b.  the  measurement 
closest  to  the  vail  vas  2 sa.  A linear  total 

temperature  distribution  vas  assumed  between  this  point 
and  the  measured  vail  temperature.  From  the  measured 
Mach  number  and  total  temperature,  the  static  temperature 
can  be  computed  using  the  energy  equation  in  the  form 

T « _ Te' 

* fF+T  (HI, 2) 

22.  Finally  the  velocity  profile  may  be  calculated  

employing  the  definition  of  Mach  number  M ■ u/a  ■ v/ftRT 
and  obtaining 
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This  velocity  profile  is  shown  in  Figure  10c  vhich  includes 
the  results  of  fcur  independent  runs.  A 1/7  - power 
profile  is  indicated  for  comparison. 

23.  8ince  nil  flow  properties  in  the  boundary  layer  are 
nor  known,  the  integrand*  in  the  equations  defining 
displacement,  momentum  and  energy  thickness,  (see  section 
71),  could  nos  be  cslculated  and  plotted.  The  integrals 
themselves  were  determined  using  a planimeter.  The  results 
are  given  in  Figure  10  and  Table  III, 2.  Comparing  these 
ratios  to  those  given  by  the  estimates  of  either  Tucker 
(references  18  and  35)  or  Eckert  (reference  24)  for  the 
insulated  flat  plate  case  at  the  same  M«o  , ve  find  the 
following  qualitative  effects  of  best  transfer  to  the  vail: 


TABLE  III, 2 


= 7 


Measured  Case: 
T-/To  .5.7 


Estimate  (reference  18  or 
24),  assuming  insulated  flat 
plate;  Pr  s 1,  1/7  - power 
profile  and  T¥/&»  - 10.8. 
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The  direction  of  deviation  of  the  measured  values 
from  those  derived  theoretically  is  expected  from 
the  increased  density  near  the  wall  in  the  experi- 
mental case  with  heat  transfer.  (Cooling  of  the  wall 
has  an  effect  qualitatively  somewhat  similar  to  suction.) 

24.  The  Reynolds  number  based  on  free  stream  conditions 
and  displacement  thickness  is  43,000.  On  the  other  hand 
an  estimate  of  an  equilvalent  flat  plate  Reynolds  number 
was  made  assuming  that  free  stream  conditions  of  the 
survey  station  would  have  prevailed  outside  the  boundary 
layer  of  an  equivalent  flat  plate.  Using  the  incompress- 
ible formula  for  boundary-layer  thickness  as  function  of 
Reynolds  number  or  that  obtained  by  Eckert  (reference  24, 
equation  37)  for  a compressible  case  with  unaltered 
velocity  profile  and  no  heat  transfer,  we  obtain  a flat 
plate  Reynolds  number  of  the  order  1 O'.  Conversely,  our 
survey  station  would  correspond  to  that  taken  on  a plate 
with  a turbulent  boundary  layer  of  about  1,5  m length 
and  our  free  stream  condition.  Our  small  tunnel  may  then 
be  taken  to  simulate  conditions  in  a large  tunnel  capable 
of  operating  with  a flat  plate  model  of  such  length 
successfully,  if  the  small  pressure  gradient  and  previous 
boundary-layer  history  do  not  affect  a fully  developed 
turbulent  bouncary  layer. 

25.  In  section  IV  it  will  be  seen  that  the  temperature 
in  the  steel  nozzls  blocks  urops  linearly  with  distance 
from  the  surface  after  a steady  state  of  supply  conditions 
and  cooling  system  operation  has  been  reached.  The  time 
history  of  such  a measurement  is  shown  on  Figure  18.  A 
cut  parallel  to  the  T - axis  through  this  figure  at  some 
time  produces  plots  of  temperature  distribution  like  that 
given  on  Figure  19.  Since  the  heat  conductivity  of  the 
nozzle  material  is  known,  the  rate  of  heat  transfer  per 
unit  area  at  the  boundary- layer  station  could  be  determined 
from 

q = A <±l  . 

n oly,  (iii, 4) 

On  the  other  hand  the  heat  transfer  coefficient  is  defined  by 


q- 


(III, 5) 


11 

UNCLASSIFIED 


- A 


E&aSagsssS! 


-mm 


1-1 


Jjg^ 


UNCLASSIFIED 
NAVQED  Report  2701 


26.  We  could  now  detenaine  the  Stanton  number  (some- 
times called  dimensionless  heat  transfer  coefficient  0^,^  ) 

q-l  __  J\  _ (cl’p&y) 

rvp  /U,_  r c ajl  f t — t \ 5 (in, 6) 

V \°°  '-■p\oo  <=o'-,e  ’*«/ 

if  we  knew  the  insulated  flat  plate  temperature,  Te. 

However,  since  Te  ££  TQ  600°K,  we  were  not  able 

to  measure  its  value  at  present.  We  therefore  calculated 
Stoo  in  two  ways.  First  we  take  the  known  approximate 
expression  for  the  turbulent  recovery  factor  (reference  26) 
r «s  ^fpr  • ^en  choose  Pr  = 0.73  based  upon 
thermodynamic  properties  near  the  wall  aud  obtain  TQ  = 539°E 
and  St « z 0.00075.  Secondly  by  using  Pr  = 1 (ana  thus 
Te  = T0  = 593°K) , we  obtain  St^  = 0.00061. 

27.  To  obtain  the  local  friction  coefficient  we  apply 
the  relation  (givsn  in  reference  26,  equation  20). 


N ^JPr  , 


(111,7) 


m 


-rr- 


which  may  be  transformed  into 


Cf  __  Qf  d‘ 

— 5~  — v-'  i r 


(III, 8) 


Using  (III, 8)  we  obtain  numerically  Cf  **  0.0012  for  both 
combinations  of  Pr  and  Te  discussed  above.  We  now  calculate 
the'  shear  stress  at  the  wall 


% = C* 


(III, 9) 


28.  Finally,  it  is  possible  to  use  this  value  for  the 
representation  of  the  velocity  profile  in  the  parameters 

. "U.  ^ 


U.+'—  7^=*  > »n<i  'U+=  ■jfrLf'iL.  (II 

(reference  27).  These  parasfeeters  are  based  on  wall 
density  and  kinematic  viscosity  since  this  plot  is  of 
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greatest  interest  near  the  wall.  Figure  3 shows  our 
measurement  with  the  region  based  upon  interpolated 
total  temperature  data  given  as  a dashed  line.  For 
comparison,  von  Karman's  incompressible  flow  semi- 
logarithaic  law  (reference  2?)  and  assumed  linear 
velocity  distribution  (1^—  "y*  } iu  the  laminar 
sublayer  are  shown.  Finally  a compressible  flow 
profile  calculated  from  equation  72  of  van  Driest *s 
analysis  (reference  25)  for  our  M «o  and  T*/T«©  is 
given.  Aside  from  other  assumptions  inherent  in 
van  Driest* s analysis,  the  calculated  profile  has  the 
boundary  condition  - -<JUo  at  y = & , matching 

the  measured  profile  at  this  point. 

29.  Summarizing,  it  may  be  stated  that  the  boundary- 
layer  characteristics  near  the  nozzle  exit  are  quite 
similar  to  those  encountered  in  incompressible  turbulent 
boundary  layers.  Although  the  heat  transfer  to  the  wall 
does  not  affect  the  velocity  distribution  materially,  it 
alters  the  relative  magnitudes  of  displacement  and  momentum 
thickness  with  respect  to  the  insulated  flat  plate  case. 
Finally,  one  cannot  discern  a laminar  sublayer  from  our 
data  possibly  because  we  were  unable  to  take  total 
temperature  measurements  in  that  range  which  is  roughly 
coincident  with  the  estimated  extent  of  the  incompressible 
laminar  sublayer. 


IV.  MEASUREMENT  OF  SURFACE  TEMPERATURES 
AND  TOTAL  HEAT  TRANSFER 

30.  Any  experiments  to  determine  recovery  temperatures, 
heat  transfer  data,  and  many  other  types  of  investigations 
rely  upon  the  accuracy  with  which  the  wind-tunnel  supply 
temperature , T0,  is  known.  The  correct  reference  Tc  is 
assured  if  the  supply  temperature  is  uniform  and  constant 
over  practically  the  entire  nozzle  inlet  cross-section. 

The  achievement  of  such  conditions,  Figure  12,  was  there- 
fore a prerequisite  to  the  experiments  described  in  the 
following  (reference  28). 

31.  Due  to  design  (cooling  system)  and  operations!  range 
(high  M,  high  T0)  of  the  hypersonic  tunnel,  experimental 
arrangements  to  obtain  wall  surface  temperatures  or 
overall  heat  transfer  data  have  to  be  different  from  those 
used  at  other  places  (reference  29).  There  is  some 

13 

UNCLASSIFIED 


z-^  --:.p  . ■ .: 
T •*  . . 


UNCLASSIFIED 
NAVOBD  Report  2701 


question  about  the  precision  that  can  be  obtained  if 
surface  temperatures  are  measured  with  a temperature 
gradient  in  the  nozzle  wall  induced  by  the  nozzle 
cooling  system.  However . the  results  obtained  show 
that  practically  steady-state  conditions  are  reached 
within  a few  minutes  after  starting  a tunnel  run,  and 
that  errors  due  to  conduction  and  radiation  are  probably 
of  the  same  order  of  magnitude  as  the  experimental 
scatter  of  the  data. 

32.  Measurement  of  Surface  Temperature.  With  the 
nozzle  cooled  by  a measured  rate  of  water  flow  through 
the  cooling  system  (see  section  lx),  a steady-state 
temperature  distribution  will  be  established  in  the 
nozzle  wall.  This  distribution  is  a function  of  the 
temperature,  pressure,  Mach  number,  coolant  flow,  etc. 
of  the  tunnel  operation.  Depending  on  these  variables 
and  the  amount  of  preheating  done  through  the  by-pass, 
it  will  take  some  time  until  practically  steady-state 
temperatures  are  attained.  If  there  are  no  serious 
non-uniformities  in  the  wall  material,  no  strong 
tempera ture  gradients  in  the  nozzle  in  the  z-direction, 
and  if  the  nozzle  wall  can  be  practically  treated  as  a 
slab  of  fixed  dimensions,  then  the  temperature  in  the 
nozzle  wall  will  drop  linearly  with  distance  from  the 
nozzle  surface.  By  measuring  the  temperatures  at 
various  depths,  it  should  then  be  possible  to  obtain 
surface  temperatures  by  linear  extrapolation . In  this 
fashion  one  avoids  the  insertion  of  measuring  elements 
that  may  disturb  the  flow  on  the  nozzle  surface. 

33.  Temperatures  of  the  nozzle  wall  can  be  measured 

at  any  of  13  points  along  the  center  line  of  each  nozzle 
block  with  four  thermocouples  imbedded  into  the  nozzle 
wall  at  various  depths  from  the  surface,  (see  Figure  2). 
To  make  installation  and  exchange  of  this  arrangement 
more  convenient,  the  thermocouples  were  inserted  in 
separate  plugs  (see  Figure  34).  These  plugs  were  ground 
to  fit  appropriate  holes  in  the  nozzle  wall.  With  the 
plugs  made  of  the  same  material  as  the  nozzle  wall,  it 
is  assumed  that  they  do  not  cause  a greater  disturbance 
of  the  tempera ture  field  than  does  the  presence  of  the 
thermocouples  itself.  To  reduce  conduction  errors  along 
the  thermocouple  wires  in  the  plugs,  they  are  led  for  a 
certain  length  behind  the  junction  parallel  to  tie  nozzle 
surface,  in  the  direction  cf  the  z-axis. 


/ -.i 
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34,  The  e.m.f.  output  of  the  thermocouples  is 
measured  with  Brown  teaperature  recorders  (model 
153X11F28A1)1  (Superscripts  refer  to  the  list  in 
section  VIII)  or  a K-2  type  Leeds  and  Northrun2 
potentiometes , dc pending  upon  the  teaperature 
difference  across  the  nozzle  wall  and  the  accuracy 
required.  With  the  Brown  recorders  the  teaperature 
of  the  individual  points  can  be  read  to±0.3QC 
(reference  28).  Using  the  K-2  type  potentiometer 
wall  temperatures  and/or  wall  temperature  differences 
can  be  read  consistently  to  within  0.01°C  during  actual 
operation  of  the  tunnel.* 

35.  Host  of  the  tests  described  here  were  made  with 
the  cooling  rate  adjusted  to  6.5  gallons  per  minute 
per  nozzle  block.  With  this  rate  practically  steady- 
state  readings  are  achieved  in  a few  minutes.  (A 
change  of  this  cooling  rate  just  gives  a displacement 
of  the  temperature  level.)  In  this  case  the  coding 
water  passed  through  two  passages  near  the  throat 
(Figure  2).  One  passage  at  the  nozzle  exit  was  partially 
open  and  all  the  remaining  ones  were  closed.  Measurements 
of  surface  tempera ture  at  Mg  = 7.6  given  here  as  example 
cover  the  temperature  range  340°K  S T0  £ 615°K 
(Figure  13)  and  the  pressure  range  9 ate  S p&  S 15  atm, 
(Figure  14).  These  measurements  were  taken  near  the  nozzle 
threat  and  it  can  be  seen  that  in  all  cases  the  nozzle 
surface  temperature  (distance  = 0 cm)  is  considerably 
below  the  tunnel  supply  teaperature  and  equal  to  or 
slightly  higher  than  room  teaperature.  An  increase  in 
supply  pressure  (Figure  14)  and  in  supply  teaperature 
(Figure  13)  causes  an  increase  of  the  steady-state  nozzle 
surface  temperature  due  to  the  increased  rate  of  heat 
transfer  for  a fixed  rate  of  coolant  flow.  Measurements  at 
given  p0  ar  T0  at  nine  stations  along  the  center  lino  of  the 


*This  accuracy  requires  the  use  of  selected,  calibrated 
thermocouple  wire,  same  wire  material  from  hot  junction 
to  ice  bath,  Cu-wire  from  ice  bath  to  potentiometer  and 
galvanometer  (i.e.  the  same  as  used  within  these  in- 
struments) , careful  aad  repeated  balancing  of  the 
potentiometer  circuit,  achieved  with  Leeds  and  Northrup 
standard  cell  and  a Willard  constant  potential  storage 
battery),  and  finally  a stable  reference  ice  bath. 
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nozzle  surface  (Figures  15  and  16)  were  evaluated  to 
extrapolate  nozzle  surface  temperatures,  wall  temperature 
gradients  (Figures  15  and  16)  and  to  construct  isotheraes 
through  the  nozzle  wall,  (Figure  17).  Finally  measure- 
ments of  wall  temperatures  at  the  boundary-layer  survey 
ctiiticn  (see  Figure  2)  were  repeated  with  higher  precision 
using  the  K-2  type  potentiometer.  The  results  indicate  a 
slow  increase  of  the  temperature  at  each  thermocouple 
location  with  time,  (Figure  18).  This  is  possibly  due 
to  the  fact  that  ei«act  steady-state  conditions  cannot  be 
reached,  as  long  as  radiative  heat  transfer  from  the  un- 
cooled sidewalls  whose  temperature  increases  slowly 
during  the  test  is  present.  However,  the  fact  that  this 
temperature  change  is  only  about  0.7°K  in  3/4  hours  permits 
us  to  consider  operation  as  practically  stationary. 

36  Assuming  no  longitudinal  heat  conduction  in  the 
nozzle  wall,  complete  insulation  of  the  nozzle  block  from 
sidewall  and  test  section  frame  (a  fact  that  is  very 
nearly  true),  absence  of  radiation  effects,  etc., 
temperature-depth  curves  such  as  shown  may  be  used  to 
determine  the  local  rate  of  heat  transfer.  (Tbe  rate 
of  local  heat  transfer  per  unit  area  is  given  by  the 
product  of  heat  conductivity  of  the  material  times  the 
temperature  gradient)  (see  equation  III ,4).  Figure  19  is 
a cross-plot  of  Figure  18  averaged  for  10  different 
abscissa  values  to  eliminate  the  small  time  effect.  With 
k3*s£i  - 50.4  keal/h  si  °K  and  at/a»  = 48.1  °S/m  from 

Figure  19,  we  find  a local  rate  of  heat  transfer  per 
unit  area  q « 2424 ±73  keal/m2  hr  at  the  boundary-layer 
survey  station  ( Moo  « 7) . This  rate  of  heat  transfer 
at  the  test  section  is  very  low  in  contrast  to  the  throat 
region,  mainly  due  to  the  low  test  section  density. 

(Compare  e.g.  the  heat  output  of  675  keal/m2  h from  a 
steam  radiator  under  standard  conditions,  i.e.  room 
temperature  24°C,  steam  temperature  110°C,  reference 
30,  p.  91). 

37.  Tests  to  determine  tbe  total  rate  of  heat  transfer 
require  again  a constant  pQ,  a constant  and  uniform  T0 
and  u oteady-atate  temperature  distribution  across  the 
nozzle  walls.  The  quantities  to  be  measured  are  the 
coolant  rate , the  temperature  of  the  incoming  or  outgoing 
coolant,  and  the  teaperature  increase  of  the  coolant. 

The  cooling  rate  is  measured  with  a Flowraior,  (see  section 
IX),  after  the  coolant  pressure  is  equalized  in  the  lines 
to  both  nozzle  blocks.  The  temperature  of  the  incoming 
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Nr? 


COOLING  RATE 
(GALLONS/MIN) 


WATER 

TEMPERATURE 

INCREASE 

(°K) 


HEAT  FLOW  RATE 
(6.27xlO~2x£TxCOOLING  RATE) 
(103  KCAL/HR) 


Ms  = 8.C  PQ  = 20  ATM  T^ 600°K 


6.6  ± 0.8% 
6.5  ± 0.9% 

5.3  ± 0.9% 
5.2  ±1.0% 
3 . 9 ± 0.8% 
3.8  ± 0 9% 
2.6±  1.9% 
2.5 ± 0.8% 

1.4  ± 1.5% 


1 . 50  ± 0.5% 
1.50  ± 0.5% 
1.84  ± 0.6% 
1 .89  ± 0.8% 
2.47±  0.8% 


2.55±  0.7% 


3.50±  0.6% 
3.68  ± 0.9% 
6.30  ± 0.6% 


2.25±  0.054 
2.23  ± 0.054 
2.20 ± 0.036 
2.22  dr  0.040 
2.17±  0.036 
2.19  ± 0.036 
2.07  ± 0.050 
2.04  ± 0.036 
1.92  i:  0.040 


Mg  = 7.6 

pQ  = 10  ATM 

Tq  = 593°K 

6.3  ± 2.3% 

0.92  ± 2.2% 

1.32±  0.06 

It-  r 


TABLE  I -IV  SUMMARY  OF  RESULTS  OF  TOTAL  HEAT 
TRANSFER  MEASUREMENTS 
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and  outgoing  water  is  measured  in  the  center  of  the 
inlet  duct  and  in  the  center  of  both  outlet  ducts, 
respectively,  with  30  gauge  I.C.  thermocouples.  The 
thermocouples. are  immersed  into  the  water  for  about 
2.5  cm,  with  the  unsupported  portions  coated  with 
Glyptai  and  scaled  into  the  support.  Because  of  the 
range  of  temperature  differences  involved  (order  of 
1°  to  10°C),  the  thermocouple  e.m.f.  is  measured  with 
the  K-2  type  potentiometer.  Measurements  at  various 
cooling  rates,  1.5  gallons/min  to  6.5  gallons/rain 
per  nozzle  block,  were  made  for  M«  = 8.0  and  pc  = 20  atm, 

Tc  z 600°K  (Figure  20)  the  conditions  for  which  the 
heat  transfer  will  be  calculated  in  the  next  chapter: 

The  results  are  given  in  Figure  21  and  Table  I. 

38.  The  errors  given  in  Table  I are  those  resulting 
from  the  scatter  of  the  experimental  data.  Estimates 
show  that  in  comparison  to  them  errors  due  to  conductive 
and  radiative  heat  transfer  to  and  from  the  nozzle  surface 
and  surrounding  boundaries  can  probably  be  neglected. 

The  nozzle  blocks  are  well  insulated  from  most  of  the 
surrounding  metal  parts  by  an  air  space,  silicone  rubber 
seals,  and  only  a small  band  of  metal-to-metal  contact. 

(The  temperature  gradient  across  the  metal-to-net&l 
contact  is  comparable  to  that  existing  in  the  subsonic 
part  of  the  nozzle,  see  Figure  17).  Conduction  errors 
are  therefore  believed  to  be  negligibly  small.  Radiative 
heat  transfer  from  the  ’-needed  tunnel  walls  can  also  be 
neglected  if  the  walls  are  kept  smooth  and  ueoxidissd. 

For  an  order  of  magnitude  estimate  of  the  error  we  assume 
a pessimistic  case  where  the  sidewalls  of  the  tunnel  over 
their  entire  lengths  are  at  supply  air  temperature  while 
the  nozzle  walls  are  kept  at  room  temperature.  Ne  assume 
further  that  the  total  energy  radiated  from  the  sidewalls 
is  transferred  to  the  nozzle  neglecting  the  geometry. 

Using  Qrad  = A cr  £ (Tj  - T24)  (reference  31,  p.  215), 

with  A - 360  cm*2  for  the  sidewall  area,  th8  Stefan-Boltzmann 
constant  <r  m 1.37  x 10”-*-2  cal/sec  cm2  0^4^  Tsidewall  = 

600°K,  Tnozzle  = 300°S,  we  obtain  Qr&d  °f  22  keal/hr  if 
the  emittance  of  the  sidewalls  is  0.1.  This  upper  bound 
of  the  radiative  heat  transfer  from  the  sidewall  to  the 
nozzle  surface  amounts  to  only  about  1 per  cent  of  the 
total  measured  heat  transfer  across  the  nozzle  wail, 

(Table  I),  and  say  therefore  be  neglected. 
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V.  THEORETICAL  HEAT  TRANSFER  CALCULATIONS  AND 
COMPARISON  WITH  MEASUREMENTS 

39.  Since  the  air  flowing  through  the  tunnel  is  heated 
and  the  nozzle  walls  are  cooled,  there  is  a continual 
flow  of  heat  from  the  air , through  the  nozzle  wall , 
into  the  cooling  water.  A steady-state  heat  transfer 
analysis  can  be  applied  to  this  problem  (see  section  IV) . 
This  analysis  can  be  subdivided  into  the  following  parts. 

(a)  The  heat  transfer  by  forced  convection  from 
a high  velocity,  compressible  fluid  to  a 
plane  wall. 

(b)  The  flow  of  heat  by  conduction  in  a homogeneous 
solid  of  a particular  geometry. 

(c)  Forced  convection  to  water  flowing  through  tubes. 

(d)  The  transfer  of  heat  by  radiation  from  the 
uncooled  sidewalls  to  the  nozzle  walls. 

The  amount  of  heat  transferred  by  (d)  may  be  neglected 
(see  section  IV) . 

A.  Overall  Heat  Transfer 

40.  If  we  at  suae  that  there  are  effective  adiabatic 
planes  midway  between  the  transverse  cooling  passages 
(Figure  22)  and  that  all  the  heat  which  enters  the  nozzle 
wall  in  the  length  A L is  removed  by  the  cooling  water, 
the  overall  heat  transfer  can  be  defined  by  the  equation 


X) 


(v,i) 


Furthermore , we  have 


_i_  + + _L_ 

-A.  Ac  Aar  A„ 


(V,2) 


from,  for  example,  reference  32,  page  136.  The 
quantities  on  the  right-hand  side  of  equation  (V,l) 
represent  respectively  the  resistance  to  heat  flow  of 
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the  cooling  water,  the  nozzle  vail,  and  the  air  streaa; 
they  are  discussed  in  detail  below. 


41.  Heat  transfer  to  cooling  water.  The  heat  transferred 
to  the  cooling  water  per  unit  area  of  nozzle  surface  is 
given  by 


L.- 


1 

Aw  rv#- 

la  calculating  bc  we  assume 


(7,3) 


(a)  that  the  entire  surface  of  the  cooling 
passage  is  At  the  saae  temperature  T^, 

(b)  that  the  teaperature  rise  of  the  water 
is  snail  enough  to  neglect  the  change 
of  u,  i,  and  Pr  with  teaperature,  and 

(c)  that  the  water  flow  is  turbulent  (true 
for  our  Reynolds  nuabers).  Ve  select 
the  heat-transfer  correlation  based  upon 
experiments!  results  for  turbulent  flow 
through  smooth  pipes  suggested  in 
reference  32,  page  167. 


°*8  r>  0 «+ 


Ki.  s=.  n.oifp  p 

! i u,  v <-  i i\  Xy" 


(V,4) 


Substituting  the  following  values  of  fluid  properties 
taken  at  14®C  (the  teaperature  of  our  cooling  water) 


Ms  1.J.5  x 10“^  kg/sec  a, 
k s IwO  x 10"4  kcal/sec 
D « .007S2  a 


a ° K,  and 


in  equation  (?,4)  gives 

c ~ 7-^0  Wr 


(7,5) 


It  should  be  noted  that  in  equation  (7,5)  W refers  to 
the  weight  of  water  entering  the  center  of  the  nozzle 
block  (see  section  X-X,  Figure  3).  Saif  of  this  water 
flows  through  each  side  of  the  cooling  passage. 
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Substituting  equation  (V,5)  into  equation  (V,3)  with 
A L = 0.0190  a gives 


~A~  = /0'°  ^°’(Tb- Z) 


(V,S) 


42.  Heat  conduction  in  nozzle  wall.  The  heat  transferred 
through  a homogeneous  wall  of  constant  thickness  is  given 
by 


r 


. A A&T 


(y,7) 


By  analogy  the  heat  transferred  through  the  nozzle  vaJl 
can  be  formally  defined  as 


(V,7a) 


where  the  effects  of  the  cool&nt  tube  size  and  position 
on  the  heat  flow  pattern  in  the  wall  are  put  into  an 
"effective  conductivity"  ke.  Assisting  the  cooling  passage 
to  act  as  a heat  sink  and  using  the  method  of  images, 


_ _ fv  m 

_ // Ji-  ~~z- 1 AL  I f 

Equation  (Y,8)  is  shown  graphically  in  Figure  23.  In 
our  nozzle  b = 0.0159  m and  from  Figure  23,  k/fc0  = 0.915. 
It  is  sufficiently  accurate  to  use  a a®an  value  of  k. 
Selecting  440°K  as  a mean  vail  temperature  gives  •, 

k = 0.012  kcal/sec  m °K.  and  substituting  these  values 
into  equation  (V,7a)  gives 


*0‘82S(Tl~mmTb)  cvr 

•n^ 

4S.  Heat  transfer  from  air  tc  nozzle  wall.  The  heat 
transferred  to” the  nozzIcTwall  is  defined as 


t-JJl-X) 

The  calculation  of  ha  is  a major  problem  in  itself  and 
will  be  considered  in  detail  later  on. 
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44.  Equations  (V,6),  (V,9),  and  (V,10)  giving 
alternative  expressions  for  q can  be  considered  a 
system  of  three  equations  in  the  three  unknowns  q, 

Tv,  Tb  and  could  be  solved  analytically  for  particular 
values  of  Te,  Tc,  V and  h . Alternatively,  a graphical 
form  of  trial-and-error  solution  using  Figure  24  may  be 

useful  bauc6  av  19  ubukii y quicker  Than  an  analytical 
solution  and  gives  a better  feeling  for  the  effects  of 
changes  in  the  values  of  the  parameters.  To  use  this 
graph  a value  of  q is  assumed  and,  for  the  prescribed 
values  of  and  h , this  fixes  T~  - Tw,  Tw  - Tk,  and 
Tb  ~ Tc.  The  sum  of  these  three  temperature  differences 
must  equal  the  prescribed  T_  - Tc.  If  they  do  not,  r, 
new  value  of  q is  chosen  ana  the  procedure  repeated 
until  agreement  is  reached. 

B.  Calculation  of  ha 

45.  The  calculation  of  h&  must  be  based  upon  the  use 
of  the  boundary-layer  concept, i.e.  the  effects  of  heat 
transfer  to  the  walls  are  confined  to  the  portion  of  air 
flowing  near  the  walls  leaving  a core  of  air  in  the  center 
of  the  tunnel  having  the  pressure  and  temperature  corre- 
sponding tc  an  isentroplc  expansion  from  supply  conditions. 
Unfortunately,  it  has  not  been  possible  to  experimentally 
determine  the  point  of  transition  from  a laminar  to  a 
turbulent  boundary  layer  in  our  wedge  nozzle.  It  is 
known  that  the  boundary  layer  at  the  exit  of  the  nozzle 

is  turbulent  (see  section  III);  however,  the  strong 
favorable  pressure  gradient  and  the  extraction  of  heat 
from  the  boundary  layer  are  likely  to  delay  transition 
(reference  33) , and  it  is  quite  possible  that  the  boundary 
layer  is  laminar  at  the  throat.  For  these  reasons, 
calculations  will  be  made  assuming  either  a laminar  tor  a 
turbulent  boundary  layer  throughout  for  a typical  case 
M8  - 8.0,  p0  = 20  atm.,  and  TQ  * 600°K. 

43.  The  velocity  distribution  along  the  edge  of  the 
boundary  layer  is  given  with  sufficient  accuracy  by 
assuming  one-dimensional  flow  based  upon  the  geometric 
area  (see  section  II) . This  velocity  distribution  is 
shown  in  Figure  25  together  with  an  outline  of  the  nozzle 
walls.  It  can  be  seen  that  most  of  the  velocity  change 
takes  place  near  the  throat.  This  is  shown  again  by 
plotting  the  local  value  of  air  acceleration  along  the„ 
nozzle  (Figure  26);  the  peak  value  is  approximately  10'  g's . 
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47.  The  simplest  approach  to  the  boundary- layer  problem 
is  to  neglect  the  pressure  gradient  (and  any  temperature 
gradient)  along  the  nozzle.  The  nozzle  is  then  hypo- 
thetically divided  into  sections,  each  of  which  has  a 
different  but  constant  velocity  to  which  the  flat  plate 
(i.e,,  zero  pressure  gradient)  boundary- layer  laws  are 
applied.  This  was  done  by  calculating  local  values  of 
Reynolds  number  along  the  nozzle.  The  length  used  in 
computing  the  Reynolds  number  was  measured  from  the 
nozzle  inlet  screen  (the  point  shown  as  0 meters  in 
Figure  25) . The  local  skin  friction  coefficients  were 
found  from 

Of  = 0.664  Itf1/2  (laminar) 

or  1/5 

Cf  0.059  Re1'  (turbulent) 

depending  upon  the  type  of  flow  assumed.  The  heat 
transfer  was  then  assumed  to  be  given  by  the  modified 
Reynolds*  analogy  (see  paragraph  27) 


(V,ll> 


where  St  is  based  upon  free-stream  conditions.  The 
heat  transfer  coefficient  can  now  be  found  from 


A = St  c* A*.  = st  -*?{*%-  -4- 

~ A 


(V, 12) 


Using  equation  (V,ll),  c_  s 0.247  kcal/kg  °C  and  Pr  = 0.68 
(for  t*  s 500°I),  " 

a t 

\ Cf  X (V,13) 

is  obtained.  Since  A*/A  varies  much  more  rapidly  with 
distance  along  the  nozzle  theic^,  the  heat  transfer 
coefficient  will  have  its  maximum  value  near  the  nozzle 
throat  (since  A* /A  is  greatest  at  M = 1) . This  is  shown 
in  the  graphs  of  ha  vs.  x,  assuming  either  turbulent  or 
laminar  flow,  given  in  Figures  27  and  28.  Since  the  value 
of  ha  given  by  equation  (V,13)  approaches  infinity  as  x 
approaches  0,  it  was  felt  that  as  our  initial  point  was 
chosen  somewhat  arbitrarily  it  would  be  undesirable  to  use 
equation  (V,13)  too  near  it.  Instead  the  value  of  ha  from 
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x : 0 to  r.  : 0.01S  m (Ji  = 0.01)  ©as  arbitrarily  taken 
as  the  value  at  the  latter  point.  This  is  also  done 
on  all  subsequent  graphs  of  ha  vs.  x. 

48.  Because  the  analysis  given  above  does  not  consider 
the  effects  of  long! luuinai  pressure  gradients  on  the 
growth  of  the  boundary  layer  and  on  ha,  its  validity 
near  the  nozzle  throat  (where  there  is  a strong  pressure 
gradient)  is  doubtful.  For  this  reason  further  cal- 
culations were  made  by  methods  which  include  some  of 
the  effects  of  pressure  gradients.  Again  separate 
calculations  were  made  assuming  either  a completely 
laminar  or  completely  turbulent  boundary  layer. 


49.  In  both  cases  the  methods  used  were  based  upon 
the  momentum  - integral  equation  (reference  34) . 


«/9  , f 


l JM 
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’-f 


/ 

& 


dot 


6* 


(v, 14) 


For  the  laminar  boundary  layer  the  method  of  Howarth 
(reference  35)  as  modified  by  Low  (reference  36)  was 
used.  This  method  is  based  upon  the  assumption  of  a 
one  parameter  family  of  boundary-layer  velocity  profiles 
using  the  Pohlhausen  parameter  ^ which  is  a function 
of  the  local  boundary  layer  thickness  and  pressure 
0-adisnt.  Since  in  Our  case,  values  of  larger  than 

the  maximum  given  in  reference  36  were  obtained,  the 
method  was  extended  using  some  of  the  results  given  by 
Thwaites  (reference  37).  Low's  analysis  also  includes 
the  following  assumptions:  (a)  constant  Ti  across  the 

boundary  layer  (the  result  of  assuming  zero  heat  transfer 

Tnd«?«i,^  = C (T/T»  ),  and  (c)  constant 

T_  along  the  nozzle  surface  (implied  by  assumption  (a)). 
Since  this  method  involved  a changing  velocity  profile 
both  the  growth  of  6 with  x and  the  value  of  for  a 
gxven  6 are  affected  by  the  value  of  the  local  pressure 
gradient.  Equations  (V,ll)  and  (V,12)  ©ere  used  to 
compute  hft  since  no  better  relation  exists  between  heat 
transfer  and  skin  friction  for  compressible  flow.  The 
variation  of  0 and  ha  with  x assuming  laminar  flow  is 
shown  in  Figure**.  29  and  30. 


50.  For  the  assumption  of  turbulent  flow,  the  method 
of  Tucker  (reference  38)  was  used,  in  this  analysis  a 


i 

i 
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constant  boundary-layer  profile  given  by  the  power 
law  i 

sl  ~ (y\  n 

U; 


(V , 17  ) 


4 n «<<■  A/t  MM/)  4L  A mI»4  m JmJ  A A J *-  ■*  — - SSSnniMM  u M 

AM  **MV « } MMM  VMV  OMAU  AA  AV  W AS/U  AO  tUSDUiUUU  AS 


Tv  __  P0i31 

fUa  ~~  Re 1/7  <Y»18) 

The  length  Reynolds  number  along  the  nozzle  surface  is 
of  the  order  o£  10® , and  consequently  a value  of , N a 7 
was  used.  In  equation  (V.18),  f and  Re  were  based 
upon  stagnation  temperature.  Since  reference  38  assumes 
zero  heat  transfer  and  i "turbulent  Prandtl  number ” equal 
to  1,  T£  is  constant  across  the  boundary  layer  and  equal 
to  T«.  The  tabulated  values  given  in  reference  38  begin 
at  M~e  0.1.  For  the  initial  growth  of  the  boundary  layer 
where  M<0.1  an  incompressible  form  of  the  momentum 
integral  equation  was  used.  Assuming 


Tw*  _ 

Oi 


rip  i i Q Tw*  Q'03.53. 

fflfc.o,  H-f  «< 

and  integrating  equation  (?,!■£>  gives 

.¥*  , 3ae,  , /..  x*-''?!. 

» = v--1  vr-  y u»  dx + &0t  (v,i9) 


from  which  the  growth  of  the  boundary  layer  can  be 
computed  (the  subscript  0 referes  to  conditions  at  the 
chosen  initial  point) . Since  0O  appears  only  in  the 
last  term  of  equation  (V,19),  it  is  evident  that  since 
u«o(/u«o  decreases  rapidly  downstream  of  the  starting 
point  (see  Figure  25) , the  effect  of  the  initial  value 
of  -9  becomes  negligible.  This  means  that  the  calculated 
thickness  of  the  boundary  layer  at  the  nozzle  throat 
is  practically  independent  of  the  assumed  distance 
from  the  throat  to  the  boundary-layer  starting  point.  The 
values  of  fca  were  again  computed  using  equations  (Y,ll) 
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and  (V,12).  The  variation  of  0 and  ha  with  x assuming 
turbulent  flow  is  shown  in  Figures  31  and  32. 

51 . Poundary-layer  thickness  at  nozzle  throat.  It 
should  be  note#  that  'despite  our  small  nozzle  throat 
opening  (equal  to  0.025*’  at  Mg  - S . 0)  the  boundary 
layer  at  this  point  is  still  relatively  thin.  Assuming 
turbulent  flow,  the  calculation  of  paragraph  50  gives 


2 x S*  (throat) 
throat  opening 


2 x 0.00016  * 0.013 
0.025 


52.  All  the  heat-transfer  coefficient  distributions 
which  have  been  calculated  are  shown  in  Figure  23 
where  they  can  be  more  easily  compared.  The  average 
heat-transfer  coefficient  E.  u£S  bovii  ck  iuuifc^60  lor  6&CD 
case,  and  the  values  are  tabulated  in  Table  V-I  together 
with  the  maximm  value  of  h*  (in  units  of  kcal/sec  a2  °K). 
The  value  of  ^ is  important  in  determining  the  overall 


TABLE  V-I 


Case 

*a 

hfc  (max.) 

Laminar  with  pressure  gradient 

0.033 

2.0 

Laminar  without  pressure  gradient 

. 0.0091 

0.18 

Turbulent  with  pressure  gradient 

0.067 

3.3 

Turbulent  without  pressure  gradient 

0.060 

1.5 

cooling  requirements  for  the  nozzle,  while  the  maximum 
value  of  lTa  is  important  in  determining  the  maximum 
temperature  on  the  nozzle  surface.  In  some  cases  the 
maximum  surface  temperature  (which  occurs  very  near  the 
throat)  say  be  more  of  a design  limitation  than  the  overall 
cooling  load.  It  can  be  seen  from  Table  V-I  that  including 
the  effect  of  pressure  gradients  in  the  laminar  boundary- 
layer  case  triples  £&  and  multiplies  ha  (max.)  by  ten. 
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Id  the  turbulent  boundary-layer  case,  including  the 
effect  of  pressure  gradients  doubles  ha  (max.)  but 
only  increases  Fa  slightly. 

C.  Estimation  of  at  Nozzle  Throat 

53.  The  nozzle  surface  temperature  couia  ts 
calculated  using  Figure  24  *nd  the  values  of  ha  found 
in  the  preceding  section.  Near  the  throat,  however,  hn 
varies  rapidly  while  Figure  24  is  based  upon  having  a 
constant  value  of  h3  over  the  length  AL.  Nevertheless, 
conservative  estimates  of  Tw  at  the  throat  (in  the  sense 
of  giving  too  high  a value  of  Tw)  have  been  made  using 
the  four  values  of  ha  (max.)  previously  calculated.  The 
calculation  is  also  conservative  because  it  neglects  the 
conduction  of  heat  in  the  nozzle  wall  parallel  to  the 
surface  away  from  the  maximum  temperature  point.  The 
calculation  was  made  for  = 587^K  (based  upon  a recovery 
factor  of  0.88),  Tc  = 287°K,  and  Wc  = 3 eal/mir.  (Lhe 
estimated  coolant  flow  rate  through  the  cooling  passages 
on  each  side  of  the  throat;  see  section  IV).  The  computed 
values  of  Tw  and  Tb  at  the  nozzle  throat  are  shown  in 
Table  V-II. 


TABLE  V-II 


Except  for  the  assumption  of  laminar  flow  without 
pressure  gradient.  Table  V-II  shows  that  most  of  the 
temperature  drop  takes  place  in  the  steel  nozzle  block. 
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D.  Comparison  of  Calculated  Values  ox  Heat  Transfer 
with  Measurements. 

54.  Using  the  measured  nozzle  surface  temperature 
(see  section  IV)  and  the  estimated  values  of  T™  st 
the  threat,  the  rate  heat  transfer  to  the  nozzle 
blech  can  oe  calculated  from 


a =-4  A (Te.-t) 


(V ,20) 


where  Tg  - Tw  is  an  average  temperature  difference 
along  the  nozzle  surface.  For  a coolant  flow  rate  of 
6,5  gal/min  a value  of  Te  - *TW  equal  to  225°K  was  used 
(where  Ts  is  based  upon  a recovery  factor  of  0.88). 

The  values  of  Q calculated  from  equation  (V,20)  and 
the  measured  value  of  Q (Figure  21)  to  one  nozzle  block 
for  a coolant  rate  of  6.5  gal/min  per  nozzle  block  are 
shown  in  Table  V-III. 


TABLE  V-III 


Q (kcal/sec) 


Laminar  with  pressure  gradient 
Laminar  without  pressure  gradient 
Turbulent  with  pressure  gradient 
Turbulent  without  pressure  gradient 


Experimental  value 


it  can  be  seen  that  the  calculated  value  of  Q based 
upon  the  assumption  of  a completely  laminar  boundary- 
layer  and  including  the  effects  of  pressure  gradients 
is  closest  to  and  lower  than  the  measured  value.  This 
is  the  condition  which  should  prevail  if  transition 
from  a laminar  to  a turbulent  boundary  layer  occurred 
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several  centimeters  downstream  of  the  nozzle  throat 
since  the  difference  in  heat  transfer  between  laminar 
and  turbulent  flow  downstream  of  such  a transition 
point  would  be  small  compared  to  the  heat  transfer 
from  the  laminar  boundary-layer  upstream  of  the  point. 

E.  Estimation  of  Heat  Transfer  Rate  from  Boundary- 
Layer  Measurements . 

55.  It  is  interesting  to  note  that  an  estimate  of 
the  total  heat  transfer  rate  to  the  nozzle  walls  for 
lls  = 7.6,  p0  s 10  atm.,  and  TQ  = 593°k  can  be  obtained 
from  the  boundary-layer  measurements  given  in  section 
III.  This  heat  transfer  rate  can  then  be  compared 
with  the  value  (given  in  section  IV)  obtained  from  the 
cooling  water  temperature  rise  for  the  same  conditions. 
Consider  the  growth  of  a boundary  layer  from  stations 
1 to  2 as  shown  in  the  sketch  below. 


f V ' 

i 


.ToeAMU^L 

& 


It  is  assumed  that  the  streamline  aE  is  always  outside 
of  the  velocity  and  temperature  boundary  layers.  It  is 
also  assumed  that  is  constant  and  equal  to  T0  across 
the  boundary  layer  at  station  1,  and  that  cp  is  constant. 
From  the  conservation  of  energy  law,  the  he&t  transfer  Q 
to  the  wall  between  stations  1 and  2 must  equal  the 
difference  in  the  energy  flowing  p^st  station  1 and 
station  2, 

<3  * v Tojr,  * Jy-  f%  “i(Tl  J y (v.2i) 


but  from  the  continuity  law 


(V,22) 
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combining  equations  (V,21)  and  (V, 22)  gives 

. Qr-"1-'  (r  • f ' H1/!  / h: 


( A-iK  <4]  J1 


(V , 23) 


rntm  4-  4 rm  / IT  /ton  Via  ttforl  f rt  S na  on  onoy*  rvj» 

WAS/l*  \ 1 • •**//  Vt**4  W W Mt^VW  V V MWA 


AI^UU  VAVU  \»  t ***//  V*»M  ~ 

thickness  0 such  that 


(v ,24) 


where  ^ \ r 1 

I felfx -£jJy 

In  applying  equation  (V,24)  to  our  nozzle,  station  1 
was  taken  at  the  nozzle  inlet  and  station  2 at  the 
boundary-layer  measurement  station  1 inch  upstream  of 
the  nozzle  exit.  The  heat  transferred  to  the  1 inch 
portion  of  the  nozzle  block  downstream  of  this  station 
is  negligible  compared  to  the  total.  A more  serious 
error  of  undetermined  magnitude  results  from  neglecting 

AU/V  ~ 9 rot  ~ 4-U«  n An  4-  4 

l,l(V  VU«  VAVU  VA  A Q A^AUgO  V lie  UVUMAV  **»*V  O VIA  b*VM  , 

Consequently,  the  value  of  Q obtained  must  be  considered 
approximate.  A comparison  of  Q (for  one  nozzle  block) 
obtained  from  equation  (V,24)  and  the  value  of  ^ 
measured  in  section  III , and  from  the  cooling  water 
temperature  rise  (see  section  IV)  is  shown  in  table  Y-IV. 


7^' 

( C4, 

\Coa  oq 


TABLE  V-IV 


Type  of  Measurement  of  Qfkcal/sec) 

Heat  Transfer  Rate 

Boundary-layer  energy  deficiency  0.33 

Cooling  water  temperature  rise  0.36 


The  agreement  between  these  two  measurements  is  probably 
within  the  experimental  accuracy. 
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Notation 

flow  area 


A L 


surface  area  of  coolant  tube  ( =7rD_) 

w 

nozzle  surface  area  per  coolant  tube  ( =<AL  x 

perpendicular  distance  froa  nozzle  surface  to 
center  of  coolant  tube 

local  skin-friction  coefficient 

specific  beat  at  constant  pressure 

disaster  of  coolant  tube 

local  heat-transfer  coefficient  for  the  air 

heat-transfer  coefficient  fcr  the  coolsmt 

average  air  heat-transfer  coefficient 

> 

equals  %*/■& 
thermal  conductivity 


•'effective  thermal  conductivity”  defined  by 
Eq.  (V,8) 

length 

distance  along  nozzle  surface  between 
"effective  adiabatic  planes"  (See  Fig.  22) 

Kach  number 

nozzle  Mach  number  setting 
Nusseit  number  hL/k) 
static  pressure 
supply  pressure 

pitot  pressure  (in  supersonic  flew  the 
stagnation  pressure  behind  a normal  shock) 
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Prandtl  number  (»  c AW.) 

P 

rate  of  heat  transfer  per  unit  area 
rate  of  heat  transfer 

recovery  factor  (a  [t^  - Te©j  / [t©  - T^j  ) 

gas  constant 

Reynolds  number  {-  ^ u x/ytx  ) 

Stanton  number  («  h/p  u cp) 
static  temperature 

OHnvtl  **  ^AmrtAWA  /aj  — -X- x t 

cm  Wtf*  « * AXOO  DUCMI  SLAgUKVAOU 

temperature) 

local  stagnation  temperature 

temperature  indicated  by  stagnation  temperature 
probe  (uncorrected) 

1 4 K*«-4  m>  ^amiSAtia 4>uv>a  /ml  «r>  1 a«1  . 

\M*rnrW  VOAAVM  A VVVTya  J 

temperature) 

coolant  tube  surface  temperature 
coolant  temperature 
velocity 
nozzle  width 

weight  flow  rate  of  coolant 
distance  in  flow  diroction 
distance  perpendicular  to  flow  direction 
ratio  of  specific  heats 
boundary-layer  thickness 
boundary-layer  displacement  thickness 

boundary-layer  sosentus  thickness 
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IX.  APPENDIX  . 

Construction  of  Nozzle  and  Cooling  System 

1.  The  nozzle  blocks  shown  in  Figures  1 and  2 were 
machined  from  a solid  block  of  steel  purchased  under  the 
Navy  Department  Specification  46  Sij— Grade  X.  This 
material  has  the  following  chemical  composition: 

0.31%  carbon,  0.25%  Xn,  0.045%  P,  0.055%  Si,  0.35%  Cu 
and  0.25%  Ni.  The  finish  machining  of  the  surfaces 
of  the  nozzle  in  contact  with  the  flow  was  accomplished 
by  surface  grinding  parallel  to  the  centerline  of  this 
surface.  The  surfaces  were  made  flat  to  within  0.00025 
inches,  and  smooth  within  an  average  of  16  BHS  (root- 
mean-square)  micro  inches  measured  perpendicular  to  the 
nozzle  centerline  or  flow  direction  and  5 RMS  micro 
inches  parallel  to  nozzle  centerline.  Tee  inlet  or 
subsonic  portion  and  the  .375  inch  throat  radius  portion 
were  ground  perpendicular  to  the  nozzle  using  a contour- 
shaped  wheel.  During  this  operation  the  nozzle  blocks 
were  set  up  side-by-side  so  that  each  cut  would  traverse 
both  blocks;  thus  producing  identical  throat  radii. 

During  this  set-up  the  nozzle  alignment  surface  was 
ground,  thus  providing  means  for  accurate  alignment  of 
the  nozzles  with  respect  to  each  other  when  installed 
in  the  tunnel.  After  grinding,  final  polishing  of  the 
throat  was  done  by  hand  to  produce  a 5 RMS  micro  inch 
finish.  Particular  care  was  taken  so  that  there  would 
be  no  irregularity  at  the  line  of  tangency  of  the 
supersonic  plane  surface  and  the  throat  radius. 

2.  The  method  of  aligning  and  supporting  the  nozzle  is 
shown  in  Figure  1.  The  nozzle  clamping  linkage,  which 
is  actuated  by  a jcckscrew  mechanism  at  its  lower  end, 
is  geometrically  arranged  such  that  it  exerts  one  major 
force  component  on  the  nozzle  in  the  vertical  direction 
to  seat  the  nozzle  against  its  upper  alignment  surface, 
and  another  component  to  set  it  firmly  on  its  4 adjustable 
legs.  Alignment  of  the  nozzle  crosswise  of  the  case  is 
controlled  by  two  alignment  pins  which  protrude  from 

the  bottom  surface  of  the  nozzle  and  engage  in  slots  in 
the  case. 

3.  The  change  of  Mach  number  setting  is  accomplished 
by  adjusting  the  two  screw  type  legs  nearest  the  nozzle 
throat  using  feeler  gage  stock  in  either  side  of  the 
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"tv- 


nozzle  throat.  The  final  throat  is  aeasured  with 
vernier  calipers  graduated  to  one  thousandth  of  an 
inch*  The  accuracy  of  the  throat  setting  is  within 
0,0005  inches  by  checking  with  a series  of  feeler 
gages  using  a "go  and  no' go"  system.  The  accuracy 
of  the  final  opening  measurement  of  the  nozzle  is  no 
better  than  0.003  inches  because  of  the  wedge  angle 
of  the  nczzle  and  the  human  error  of  locating 
repeatably  the  final  opening. 

4.  The  distension  of  the  gasket  grooves,  as  shown 
below  were  determined  to  both  retain  the  gasket 
material  and  to  allow  the  door  to  close  to  a mstal- 
to-metal  seat  while  affecting  a seal. 


tfOZUM 

BOB?/ 


-.109  R 


-.234 

±.003 


LbiKott 

- 3\OcW 


.109  R 


\y.zi9 

pi*)  ±X3C3 

S4UGNMEW 


NOZZLE  GASSET  GROOVE 
PARALLEL  TO  FLOW 


GASKET  GROOVE  BETWEEN 
NOZZLE  BODY  & TUNNEL  CASE 


The  gasket  material  used  throughout^ ths  tunnel  is 
"Cohrlastic  Silicone  Rubber  R-8925"3  which  was  purchased 
in  1/4  inch  diameter  extruded  cord  form.  The  best 
tolerance  available  on  the  diameter  is  ± .015  inches, 
to  this  date. 

5.  A portion  of  the  tunnel  that  1b  intimately 
associated  with  the  cooling  wedge  nozzle  is  the  air 
supply  inlet  system,  (reference  28),  Figure  12.  This 
arrangement  consists  of  a baffle  system  followed  by  a 
honeycomb  in  an  insulated  chamber  for  mixing  the  air  to 
obtain  a nearly  uniform  temperature  distribution  ahead 
of  the  nozzle  threat.  A stainless  steel  screen  in  a 
"tr.onsite"4  frame  was  found  desirable  to  prevent  clogging 
of  the  throat  and  probes  by  foreign  particles,  chips  of 
the  ceramic  insulators  and  oxide  flakes,  fross  the  exposed 
elements  of  the  electric  heaters.  This  system  produces 
a nearly  uniform  supply  temperature  at  the  nozzle  throat. 
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Cooling  Sys'.jm 


1 . The  internal  cooling  passage  arrangement  of  the 
nozzle  block,  as  shown  in  Figure  3,  consists  of  a 
central  inlet  water-distribution  passage  with  34  holes 
of  0.312  inch  diameter  located  5/3  inches  from  the 
air-flow  surface.  These  holes  connect  this  inlet 
passage  to  each  of  the  two  side  outlet  passages  and 
can  be  easily  closed  (e.g.  by  removal  of  the  side 
plates  and  stopping  the  holes  with  rubber  cord),  or 
restricted  by  insertion  of  an  orifice  plate. 

2.  Tho  exterior  piping  arrangement,  as  shown  in 
Figure  1,  is  made  up  of  slandard  1 inch  copper  water 
tubes.  The  water  supply  main  is  located  beneath  the 
floor  on  the  lower  left.  It  is  controlled  by  a valve 
attached  to  a removable  "T"  handle  (not  shown) . The 
water  is  fed  to  the  right  side  of  the  tunnel  through 
the  tube  near  the  floor  on  the  far  side  of  the  tunnel. 
The  inlot  water  pressure  was  equalized  by  the  right- 
angle  valves  and  the  rate  of  flow  is  measured,  in  the 
nozzle  under  study  by  a variable  a»*ea  tjpe  riov  Meter® 
on  the  intake  pipe  shown  at  left.  The  two  exhaust 
water  pipes  coming  from  each  nozzle  converge  and  join 
immediately  outside  of  tho  tunnel . The  flow  from  each 
nozzle  is  throttled  separately  by  valves  located  at 
the  level  of  the  circular  window  on  the  far  side  of 

the  tunnel. 


3.  ’The  coolant  used  is  tap  water.  For  the  usual 
flow  rate  of  6.5  gallons  per  minute  the  inlet  water 
pressure  is  6.5  psi  and  the  outlet  pressure  is 
throttled  to  maintain  a pressure  of  1 psi. 

Thermocouple  installation 

1 ..  The  nozzle  temperatures  are  measured  by  iron- 
constantan  thermocouples  made  of  £30  gage  wire2  and 
electric  welded  without  flux  using  a Variac  and  Pyrotip 
burner6  with  a carbon  electrode.  The  thermocouples 
were  secured  in  place  using  technical  "B"  copper  cement”^ 
in  a modified  taper  pin,  shown  in  Figure  34.  The  exact 
location  of  each  of  the  four  thermocouple  junctions  in 
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each  plug  is  determined  to  0.01  millimeters  after 
assembly  by  use  of  a tool  maker’s  microscope8  (Bausch 
and  Lcmb  - Type  UK  352).  Special  care  is  taken  in 
the  manufacture  of  the  plugs  so  that  each  plug  would 
bottom  in  its  respective  hole  in  the  nozzle  block. 


Plug  # 1 

2 

3 

4 

5 

6 

Location 

from  Final 

Nozzle  . 656 

1.406 

4.406 

7.406 

9.656 

11.906 

Opening 

C ISC-65 ) 

7 8 

9 

10 

11 

12 

13 

14.156  16.4C6 

18.656 

20.156 

21.656 

23.156 

24.156 

TUNNEL  3IDES ALLS 

1.  Several  door  arrangements  are  available  for  the  tunnel. 
There  is  one  set  of  blank  steel  doors,  one  steel  door 
with  42  centerline  pressure  taps  located  one  inch  apart, 

a set  ox  doors  having  a rectaegnl  gsJ-.te  glass  window, 
and  a set  of  doors  having  3 circular  windows  (see  window 
in  Figure  1).  The  pressure  tap  door  and  the  circular 
window  doors  sure  provided  with  two  sets  of  binge  arms 
so  that  when  the  door  is  inverted  the  tap  holes  and 
windows  bisect  their  forcer  location.  The  circular 
windows  now  in  use  are  cut  from  selected  commercial 
plate  glass.  Experience  indicates,  that  while  running 
at  supply  temperatures  of  900°F  these  windows  can  be 
used  anywhere  fro?1  the  nozzle  esit  through  the  diffuser 
for  a casisus  period  of  20  minutes  without  fracture. 

A fused  quartz  window  of  similar  arrangement  is  now 
being  tested  to  increase  the  time  of  exposure. 

2.  Fhss  the  doors  are  closed  they  are  forced  by  bolting 
against  the  raised  nozzle  lips,  which  forms  the  edges  of 
flow  surface  and  against  the  metal- to-metal  bearing 
surface,  of  the  tunnel  case  (see  Figure  1 outside  of 
gasket  ) . Xu  udui  tlon  to  bolting  around  the  periphery, 
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F!G.  S3  NOZZLE  WALL-TEMPERATURES  NEAR  THROAT 
(STATION  NO.  13)  AT  VARIOUS  SUPPLY  TEMPERATURES 
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lil 


, c ~r 


• ■ ■ v m -t? 


MP 


■ , 


Wftmrn_ 


pgril 

Itagir 


NAVORD  REPORT  270! 


B$SJ 


,§p| 


ii 


i-'-jhi 


i^i 


i 


«:•  xi 


2,0.6 

ui 

O 

z 

UJ 

a: 

UJ 

tC. 4 

o 


<0.2 

UJ 

Q_ 

UJ 


at . oAfl  !i< 

Ay  ~‘°48  MM 


4 8 12 

DISTANCE  BELOW  SURFACE  (MM) 


,rl 


i&sis 


:WWi 


mm 


ERROR 


( distance;!  o.i  mm 

B 

} AT:±0.0l°K 


FIG.  19  TEMPERATURE  DROP  IN  NOZZLE  WALL 

Mg  s 7.6 , pe  ~ 21.  4 ATM  , T4>=593°K 
COOLING  RATE:  S.5  GALLQNS/MIN 
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□ NORTH  SIDE  OF  NOZZLE 
B SOUTH  SIDE  OF  NOZZLE 
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FIG.  20  TEMPERATURE  INCREASE  OF  COOLING  WATER 
AS  FUNCTION  OF  THE  COOLING  RATE  FOR 
TUNNEL  OPERATION 
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FIG.-34  INSTALLATION  OF  THERMOCOUPLE  PLUG  IN  NOZZLE  BLOCK 


